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ABSTRACT 

We have used the KPNO 4-meter Mayall telescope to image 0.83 square degrees 
of sky in six fields at high galactic latitude in six filters spanning 3000-lOOOOA to 
magnitude limits ranging from 22.1 to 23.8. We have assembled a catalog of 21,375 
stellar objects detected in the fields for use primarily in conducting a multicolor search 
for quasars. This paper describes the data reduction techniques used on the CCD data, 
the methods used to construct the stellar object catalog, and the simulations performed 
to understand its completeness and contamination. 
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1. Introduction 

Quasars, thought to be powered by the accre- 
tion of matter onto supermassive (~ 10 7 ~ 10 Mq) 
black holes at the centers of galaxies, are the 
most luminous individual objects in the universe 
(-31<M<-23). As such, they can be seen to 
greater distances than any other objects, provid- 
ing us with a unique opportunity to study the 
universe back to times when it was only ~10% of 
its current age. Quasars yield information about 
the universe through their absorption-line sys- 
tems and associated galaxies (e.g. [Steidel Sz Sar- 
gent 199^ ), through the galaxies in quasar envi- 
ronments ( |Yee fc: Ellingson 1993 , and references 
therein), and through the lensing of quasars by 
intervening systems ( [Kaiser &: Tribble 1991 ), to 



name but a few. Additionally, the mere existence 
of quasars at the observed space densities at red- 
shifts z>4, when the universe was only ~10% of 
its current age, constrains models of galaxy and 



structure formation in the universe (e.g. Turner 



199l|, |Loeb 19931, |Nusscr fc Silk 1993j [Haehnelt 
k Rces 1993| , [Katz et al. 1994|) . 



One direct way of using quasars to study the 
universe is to map the evolution of the quasar 
luminosity function (QLF) with lookback time; 
that is, the number of quasars as a function of 



absolute magnitude and redshift (e.g., Schmidt fc 
Glrecn 1983 , Hewett et al. 1993| ). Previous stud- 
ies of the QLF have shown that the space den- 
sity of the bulk of the quasar population, with - 
28<M<-23, was considerably greater in the past, 
by a factor of ~100 at z=2 relative to the cur- 
rent comoving density. This was likely due to 
the evolution of the quasar population in both 
luminosity and in density, though the exact form 



of the evolution is not well constrained ( Hewett 



it al. 199"3| ). The redshift distribution of these 



quasars may peak at z~2-3 and decrease at z>3, 
though not as swiftly as at z<2, as first proposed 
by Osmer (1982). 

Recently Warren, Hewett, & Osmer (1994; 
hereafter WHO), using a multicolor-selected sam- 



ple of 86 quasars with 2.2<z<4.5, showed that 
there is indeed a decline in the observed space 
density of quasars with z>3. This conclusion has 
been confirmed by the recent results of Schmidt, 
Schneider, & Gunn (1995), using a Lya-selected 
sample of 90 quasars with 2.75<z<4.75. Al- 
though some groups have previously found no 
evidence for such a decline (e.g. Koo & Kron 
1988 (KK88), [Giallongo fc Vagnetti 1992j) , these 
more recent studies argue that the decline is real. 
However, it should be noted that WHO cannot 
rule out that the decline they see in the observed 
quasar space density at z>3 is due to dust ob- 
scuration, as suggested by e.g. Fall & Pei (1993), 
and not reflective of a decrease in the true space 
density of quasars at high redshift. One of the 
reasons for this contention over the behavior of 
the QLF at high redshift is that each sample of 
quasars typically contains only a few such ob- 
jects. In addition, accurate determination of the 
space density of quasars at high redshift requires 
careful modelling of the selection effects present 
in each survey in order to make incompleteness 
corrections, as discussed e.g. in Warren et al. 
(1991a). 

There is obviously a need for additional carefully- 
selected samples of quasars with z>2.2. In 1990 
we initiated a deep multicolor survey for quasars 
with the intent of obtaining such a sample. Our 
survey was enabled by the arrival of the first 
2048 x 2048 CCD at the KPNO 4-m telescope; 
the detector and telescope offered at that time 
a unique combination of large areal coverage on 
the sky, spectral sensitivity across the optical 
band, and deep limiting magnitude. The sur- 
vey was conducted in six optical filters span- 
ning 3000-9000A in six different high-Galactic- 
latitude fields. The primary objective was to 
construct a well-defined stellar locus from which 
quasar candidates could be easily isolated via 
their unusual colors. However, many other in- 
teresting questions can be addressed with the 
survey. We achieved good photometric calibra- 
tion and the star-galaxy classification effects are 
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well understood, and thus the survey should be 



useful for studies of Galactic star counts (Reid 
& jVlajewski 1993j ), The large number of galax- 
ies detected in the survey also forms a valuable 
database for the study of field and emission-line 
galaxies. Our wide wavelength coverage makes 
the survey especially useful for studies of field 
galaxy color evolution and for the identification 



our classification routines to have an appear- 
ance consistent with an unresolved point source 
in a certain minimum number of images, as de- 



of rare galaxy types such as E+A galaxies (|Liu 
et gL 1996| ). 

Quasars with -25<M<-23 at z=3-5 are found 
at apparent magnitudes of roughly 20<R<23, so 
our goal was to image to R=24 (R=23.5 was 
achieved). Such a faint limiting magnitude also 
probes the upper end of the Seyfert luminosity 
function at z<3. To our limiting magnitudes, 
the survey region contains approximately 30,000 
objects, of which about half are galaxies. Based 
on previous knowledge of the apparent luminos- 
ity function of quasars ( [Majewski et al. 1993 ), it 
should contain almost 200 quasars. Our ultimate 
goal is to make available catalogs of the quasars, 
galaxies, and stars containing accurate positions, 
magnitudes, and photometric error estimates for 
all the objects. 

The survey for quasars had two principal ob- 
jectives: 1) to investigate the luminosity function 
of high-redshift quasars (z > 3) at fainter magni- 
tudes than had been done before, for which the 
absolute magnitudes were more comparable with 
previous surveys at lower redshift; and 2) to de- 
termine the luminosity function of quasars with 
z < 2.2 at faint enough magnitudes to constrain 
the debate on luminosity vs. density evolution, 
which has been unresolved for lack of data at the 
faint end of the luminosity function, where the 
differences between the two hypotheses are pro- 
nounced. 

This paper describes the CCD imaging ob- 
servations done for the survey and the data re- 
duction methods used to construct the catalog 
of stellar objects from which quasar candidates 
were chosen. (We define a stellar object as an 
object detected in our images that is found by 



tailed in §4.3 below.) Paper II ( pall et al. 1996] ) 
discusses quasar candidate selection and initial 
spectroscopic results, and makes a comparison 
with expectations from previous luminosity func- 
tion studies. Paper III ( psmer et al. 19"96 ) will 
present final spectroscopic results and a detailed 
derivation of the luminosity function. 

This paper is organized as follows. In §2 we 
describe the observations, the survey fields, and 
the filters used. We cover in §3 the main features 
of the data reduction steps. In §4 we summarize 
how objects were detected and cataloged using 
FOCAS and how the photometric measurements 
were made. We discuss in §5 the photometric 
calibration and analysis of the errors, and in §6 
the astrometric calibrations. Finally, in §7 we 
summarize the properties of the resulting catalog 
of stellar objects. 

2. Observations 



Images were obtained at the prime focus (f/2.7) 
of the KPNO 4-meter Mayall telescope with a 
Tektronix 2048x2048 CCD with 27-micron pix- 
els, giving a plate scale of O'.'529/pixel. A gain 
of 11 e-/ADU was used, yielding a full well of 
~30,000 ADU. The particular CCD we used was 
an engineering grade, thick chip with poor cos- 
metics. It had a useful area of 1860x2048 pixels 
due to an excess of variable bad columns on one 
side of the chip, and also had a large number of 
traps evident in the lower light levels of the UBV 
images. 

We used a set of six filters covering the wave- 
length range from 3000 to 9000 A. The set con- 
sisted of standard KPNO Johnson UBV filters 
(A e// =3595 A, 4364 A, and 5563 A respectively) 
plus three special red/near- IR filters: R', 175, 
and 186 (A e// =6615 A, 7425 A, and 8586 A re- 
spectively). The response curves for each filter in 
the absence of atmospheric extinction are plot- 
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ted in Figure 1. Each curve is the product of 
the quantum efficiency of the CCD, the trans- 
mission of the filter, the transmission of the 4- 
meter prime focus doublet corrector, and one re- 
flection off aluminum. The R' filter is narrower 
than the standard KPNO Harris (not Mould) R 
filter, having less of a red tail and also being less 
transmissive overall, but it has very nearly the 
same effective wavelength (6615 A) as the stan- 
dard KPNO Harris R filter (6625 A). Thus for 
most objects it is safe to treat magnitudes in this 
filter as a standard Harris R filter; only the zero- 
point of the photometric solution will differ sub- 
stantially from the solution for a Harris R filter 
(although there will be a color term for very red 
objects). Hereafter the R' filter will be referred 
to simply as R. The 175 and 186 filters each cover 
about half the wavelength range of the standard 
I filter. These three red filters were designed with 
narrower than usual bandpasses in order to im- 
prove our ability to detect quasars with redshifts 
as large as 5.5 and to distinguish z>4 quasars 
from very late-type stars via the strong excess 
from La emission. For z<5.5, the 186 filter is 
still longward of redshifted La emission, which is 
normally the strongest line in quasar spectra, and 
therefore provides a discriminant between high-z 
quasars and later-type stars. 

Successful imaging runs, affected only occa- 
sionally by light cirrus, were carried out on 1990 
August 21-24 and 1991 April 14-17, resulting in 
images being obtained in twelve CCD fields. Our 
fields are arranged on the sky in six pairs of two 
adjacent individual CCD fields (with slight over- 
laps), and were chosen to be at high galactic lat- 
itude away from known dust clouds. A list of the 
coordinates, area, and galactic extinction param- 
eters for each field is given in Table 1. 

We attempted to obtain two sky-noise-limited 
images of each field in every filter, offset by typi- 
cally ten arcseconds from each other, to help re- 
duce spurious detections of cosmic ray hits and 
bad pixels. We succeeded in obtaining at least 
two offset images except in the Ole field and the 
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Fig. 1. — The response curves of the survey in 
each filter (UBVRI75I86 from left to right). Each 
curve is the product of the quantum efficiency 
of the CCD, the transmission of the filter, the 
transmission of the 4-meter prime focus doublet 
corrector, and one reflection off aluminum. At- 
mospheric extinction has not been figured into 
these curves. 

lOe field in the B filter. Where more than two off- 
set images were obtained, the two best in terms 
of seeing, depth, and focus quality across the chip 
were selected to be the final two images used 
in the reduction process. Various parameters of 
these 125 CCD images that comprise the survey 
are listed in Table 2. 

3. Data Reduction 

The CCD images in all filters of each field 
were reduced using IRAF0. The bias level was 
subtracted from each column separately on all 
field, dark, and domeflat images. The appropri- 
ate dark current image was subtracted from each 
field image, and domeflats were used to remove 
the pixel-to-pixel sensitivity variations. A sepa- 
rate bad pixel mask needed to be made for each 
image. Each bad area was interpolated across its 
narrowest dimension. 

3 The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Optical Astronomy Observato- 
ries, which is operated by AURA, Inc., under contract to 
the National Science Foundation. 
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It was found at this point that the images still 
showed random, nonrepeatable gradients and vari- 
ations in the sky level across the chip. This 
was particularly true in the three reddest bands, 
where the variations were typically 5% peak-to- 
peak. These variations were assumed to be ran- 
dom additive variations in the sky values, and 
were removed using a two-step process. First, 
all images in a particular filter for each run were 
combined to form a "skyflat" (identical to the 
"supersky" described in Gullixson 1992) using a 
program supplied by Todd Boroson. The images 
were scaled by the mode of a central region and 
then median filtered to form the skyflat. This 
procedure was successful at eliminating objects 
completely except when fewer than 10 images 
were available to form the median. To eliminate 
artifacts from such residuals, we smoothed the 
skyflat using a variably sized boxcar filter via the 
IRAF mkskycor task. Each skyflat was then 
normalized by the mode of its central area and 
divided into the images, resulting in a reduction 
of the gradients and variations across each image 
to at most the 2% level. To eliminate these re- 
maining random variations, which could be found 
in every filter except U, we used the Faint Ob- 
ject Classification and Analysis System (focas, 
Valdcs 1982a| ) sky-fitting routine to estimate the 
local sky value for every 8x8 pixel box. The sky 
values determined by FOCAS were found to be ac- 
curate to within 1% even in areas dominated by 
saturated stars, but were offset from their true 
locations slightly due to the line-by-line process 
of sky determination used by FOCAS. To avoid 
creating troughs to one side of bright objects af- 
ter sky subtraction, we produced three additional 
FOCAS sky images, rotating the image by 90° be- 
tween each. These sky images were rotated back 
to the original orientation and medianed with the 
first sky image, discarding the highest pixel, to 
form the final sky image. Each field image had its 
own sky image subtracted from it, leaving fields 
globally flat to <0.5%. The only shortcoming 
of this method was that the outer regions of the 



largest, most extended galaxies present (typically 
20-50" in size) were sometimes confused with the 
sky, resulting in an oversubtraction of the sky 
around these galaxies. This effect is not of con- 
cern for the catalog of stellar objects. 

Next, cosmic rays were isolated and interpo- 
lated across using the IRAF task COSMICRAYS. 
The removal worked quite well, having trou- 
ble only with some many-pixel cosmic ray hits 
and almost never flagging unsaturated objects 
by mistake. To allow us to be conservative 
in our eventual selection of candidates in the 
event that the interpolation introduced errors 
into the object photometry, all cosmic ray loca- 
tions were added to the bad pixel lists for later 
cross-correlation with the object lists. 

Lastly, the images were shifted and trimmed 
so that each CCD field was on a common coor- 
dinate system to the nearest pixel (no fractional 
pixel shifts were used) . The area covered by these 
trimmed images, taking into account the addi- 
tional reduction from adjacent field overlaps and 
the 5" wide border on each image unused by FO- 
CAS, is given in Table 1. The total areal coverage 
on the sky of the survey is 2990D', or 0.83D . 

4. Detection and Cataloguing of Objects 

Because of the variations in seeing, focus, and 
pointing between the two images in each filter 
in our fields (which were often taken on different 
nights), as well as the problem of cosmic-ray con- 
tamination, we decided not to coadd the images 
to produce a single image with improved signal- 
to-noise (S/N) in each filter. Instead, we reduced 
each frame individually and improved the S/N 
by averaging the magnitudes from each frame af- 
ter they had been reduced. The main drawback 
to this approach with respect to coadding is the 
increased unreliability of object classification at 
faint magnitudes, but we are helped in this by 
having two independent classifications. We used 
FOCAS and its 'built-in' detection filter to detect 
and classify all the objects in our frames, and the 
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IRAF phot routine to find aperture magnitudes 
for them. 

A brief summary of the detection and initial 
photometry steps performed on each field follows: 

• FOCAS detections must have at least 9 contigu- 
ous pixels 2.5a above the sky, after the value of 
each pixel was convolved with the values of pix- 
els within a 5x5 grid, weighted according to the 
standard 'built-in' FOCAS detection filter 

• Objects lying on bad pixels were flagged as such 

• Objects were classified using a different em- 
pirical PSF for each of nine subsections on each 
image 

• Objects not appearing on both images in the 
filter were discarded 

• To be considered stellar in appearance, objects 
were required to be classified star in at least 
one image of half or more of the filters in which 
they were detected 

• Aperture corrections were calculated for each 
image subsection 

• Instrumental magnitudes were calculated in a 
3'.'0 diameter aperture 

• The instrumental magnitudes were calibrated 
to a standard system 

• Magnitudes from each image in each filter were 
combined, accounting for bad pixel contamina- 
tion and magnitude limits, to give a final magni- 
tude for each object. 

Below, we discuss important reduction steps 
in more detail, where necessary. 

4.1. Bad Pixels 

All objects that fell on bad pixels were flagged 
as such. 'Bad pixel' here is a broad term, in- 
cluding sites of detected cosmic ray hits, the ac- 
tual bad columns and traps determined through 
inspection, and the bad focus regions (typically 
the corners) of each image. Saturated stars were 
also considered to fall on bad pixel areas, but 
were flagged differently since the determination 
of color limits for such objects is different than 



for regular bad pixels (we measure a lower limit 
to the object's magnitude in these cases). 

4.2. Classification 

Objects were classified using the standard FO- 
CAS method and parameters, using the resolution 
classifier fully described in Valdes (1982b) and 
summarized briefly here. An empirical point- 
spread function (PSF) is constructed from iso- 
lated bright stars. From the PSF a set of clas- 
sification templates is generated, consisting basi- 
cally of a broadened or narrowed PSF. Objects 
were then classified as star, fuzzy star, galaxy, 
diffuse object, or noise spike depending on which 
template they best fitted. Due to the frequent 
variation in focus across the chip, each image was 
divided into 9 subsections in a 3x3 grid, and a dif- 
ferent empirical PSF was used for classification 
in each subsection. This method is obviously not 
as accurate as fitting a slowly varying PSF across 
the chip; but the latter method is difficult to im- 
plement with FOCAS. We included 'fuzzy stars' 
in our stellar object list to increase our ability to 
include objects where an active nucleus is seen 
along with extended light from a host, compan- 
ion, or intervening galaxy. Inspection of the lo- 
cation of objects on a plot of PSF scale vs. nu- 
clear component fraction showed that only ~10% 
of fuzzy stars have an extended component with 
amplitude significantly different from zero. This 
tells us that extended objects with strong nu- 
clear light components are relatively rare, and 
reassures us that we are not including a large 
number of galaxies with this criterion. 

It should be noted that some of the frames had 
a strong gradient in focus across them, causing 
substantially out of focus images at the corners 
of the chip. Objects in these regions were flagged 
so that the photometry from them is not used; 
fortunately, no region was out of focus in both 
offset images of the field. 
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4.3. The Stellar Object Catalog: Selec- 
tion, Completeness, and Misclassifi- 
cation 

Both offset frames in each filter were matched 
with each other, and objects not detected in both 
offset images were discarded. An object was clas- 
sified stellar in a given filter if it was classified by 
FOCAS as a star or fuzzy star in either offset frame 
of the filter. The criteria for inclusion in the final 
catalog of stellar objects were that the object be 
detected in at least three filters and classified as 
stellar in half or more of the filters in which it 
was detected. Our criteria are rather lenient in 
allowing objects classified as nonstellar in some 
images to be included in the stellar object cat- 
alog because we did not want to exclude faint 
AGN with detectable host galaxies. We need to 
understand how often nonstellar objects are mis- 
classified as stellar, and vice versa, as a function 
of magnitude in each filter in order to understand 
the contamination and completeness of our cat- 
alog. 

The ideal method of determining the com- 
pleteness and classification reliability is to add to 
the original images simulated objects of known 
morphology and a wide range of magnitude. 
These images are then processed as the original 
images were, and the completeness and classifica- 
tion reliability are straightforward to determine. 
However, the varying PSF over many of our CCD 
images makes the addition of accurately simu- 
lated objects to the data complicated. Instead, 
we opted to create entirely artificial images and 
analyze them. Since the parameters of the artifi- 
cial data can be made to match those of the real 
data very accurately, this approach should yield 
an accurate description of the real data's classifi- 
cation reliability and completeness as a function 
of magnitude. 

We used the IRAF artdata package to sim- 
ulate two images of a CCD field in each filter. 
We attempted to match the seeing, background, 
and noise level observed in the Olw field. For the 



galaxies, to subjectively match the appearance of 
the real data, a characteristic radius of 5.3" and 
a mixture of 40% de Vaucouleurs profile 'bulges' 
and 60% exponential profile 'disks' was chosen. 
Galaxies were given a power-law magnitude dis- 
tribution with slope 0.35 (i.e. n(m)~ io°- 35 n 
which matches well the slope of faint magnitude 
galaxy counts. To increase the number of faint 
stars available for studying completeness, stars 
were given a power-law magnitude distribution 
with slope 0.6, steeper than the observed dis- 
tribution. However, since this distribution does 
not generate as many bright stars as the real im- 
ages contain, it was necessary to add a grid of 
bright stars to the images to ensure that the PSF 
could be well determined in each image subsec- 
tion. Since whether or not an object is included 
in our catalog of stellar objects depends on its 
classification and magnitude in several filters, ob- 
jects of different colors had to be simulated. We 
divided the locus of stellar objects in our real 
data into five bins, each with different broad- 
band colors. The relative numbers of stars and 
galaxies in each bin was determined from the real 
data, and corresponding percentages of the simu- 
lated stars and galaxies were assigned the colors 
of each bin. Finally, following the above parame- 
ters, 3000 stars and 2500 galaxies with R= 16-24 
were randomly distributed across the field. The 
appearance and statistics of the images closely 
matched those of the real data. One minor dif- 
ference is that the PSF in the real data is not 
a perfect gaussian. However, the difference be- 
tween the real PSF and the simulated one is only 
apparent in the wings of bright stellar objects, 
and so this difference should not affect any of 
our results on faint objects. 

Our detection and cataloging programs were 
then run on the artificial images in exactly the 
same manner as on the real data, producing an 
output FOCAS catalog of stellar objects. The FO- 
CAS task artcat was used to create FOCAS cat- 
alogs of the stars and galaxies input to the ART- 
DATA task, for comparison with the output cata- 
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log. 

To determine its completeness, the output cat- 
alog of stellar objects was compared to the input 
catalog of stars. As is typical for CCD data (|Stct- 
1990| ), the completeness fraction is 100% at 



sor 



bright magnitudes, begins a slow linear decline 
at some faint threshold magnitude as the effects 
of crowding appear, and then switches to a rapid 
linear decline as the detection threshold comes to 
dominate over crowding effects. We find that this 
transition occurs at about the 90% completeness 
level, which closely approximates the 5a detec- 
tion limit (see §4.5). Similarly, the 50% com- 
pleteness limit closely approximates the 3<r de- 
tection limit. These limits are listed in Table 3, 
along with the threshold magnitude above which 
the completeness is essentially 100%. 

To determine the classification reliability, the 
output catalog of stellar objects was compared to 
the input catalog of galaxies to determine which 
supposedly stellar objects were actually misclas- 
sified galaxies. In principle, false detections due 
to noise could also show up in the catalog, but 
by requiring detection in both offset frames of at 
least three filters we completely eliminate such 
detections. 376 galaxies were misclassified as 
stars in our simulations, virtually none of them 
brighter than the completeness threshold mag- 
nitude, which we thus adopt as the threshhold 
magnitude for classification as well. To apply 
these results to the real dataset requires adjust- 
ing for the expected number of faint galaxies in 
the survey fields versus the number used in the 
simulated data, since the number of misclassi- 
fied galaxies in the catalog obviously depends on 
the number of galaxies that exist in our fields. 
We used the compilation of McLeod &: Rieke 
(1995) to determine the number of faint galax- 
ies expected in our fields clS Si function of mag- 
nitude. The slope of the N(m) relation for our 
simulated galaxies is 0.35, close enough to the 
value observed for faint optical counts that we 
only need to adjust the normalization. To ap- 
proximate the counts from R=21. 5-23.5 we used 




Fig. 2. — The upper curve is the observed cumu- 
lative number-magnitude (R band) counts of all 
21375 objects in the stellar object catalog. The 
lower curve is the expected cumulative number- 
magnitude (R band) distribution of misclassified 
galaxies in the catalog, based on scaling the re- 
sults of the artificial data simulations described 
in the text to the expected galaxy number den- 
sity and total area of the survey. At faint magni- 
tudes (R^21.5) the stellar object catalog shows 
an excess above the power-law distribution seen 
at intermediate magnitudes, which may reason- 
ably be attributed to misclassified galaxies. 

the counts from Gunn r=22-24 and assumed that 
those counts are 100% complete. In our 0.0713D 
simulated field, the input catalog contained 1761 
galaxies from R=21.5-23.5, whereas ~50% more 
(2797) are expected in an equal area of real data. 
Thus to make our simulated field accurately re- 
flect the number density of faint galaxies, we ad- 
just the total number of contaminating galaxies 
upwards by 50%, from 376 to 564. In our en- 
tire survey area of 0.83D , we then expect 6565 
spurious objects in the stellar object catalog. 

Using this scaling and the magnitude distri- 
bution of the misclassified galaxies in our simu- 
lation, in Figure 2 we compare (in the R band) 
the expected number of contaminating galaxies 
to the total number-magnitude counts of the stel- 
lar object catalog. At faint magnitudes (R;>21.5) 
the stellar object catalog shows an excess above 
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the power-law distribution seen at intermedi- 
ate magnitudes. This may reasonably be at- 
tributed to misclassified galaxies. A further 30% 
increase in the number of contaminating galax- 
ies, to 8535, is required to provide a good match 
to the number-magnitude distribution of this ex- 
cess. 

What might have caused our simulations to 
under predict the number of misclassified galax- 
ies? There are some effects we have neglected 
in the simulations, but most of them should not 
alter our conclusions significantly: 

• We have not included cosmic rays in the ar- 
tificial data, although this capability does exist 
within artdata. Cosmic rays would effectively 
act as another source of noise, and since our se- 
lection criteria are very efficient at keeping false 
detections due to noise out of the catalog, their 
effect is negligible. 

• We have not included the effects of interpolat- 



in g! over bad pixels, although this is not difficult 



to do either. Interpolation might cause objects 
to fall below the detection threshold or be mis- 
classified. Since the area affected by bad pixels is 
very small and its location varies from image to 
image, and since our criteria for inclusion in the 
final catalog uses information from many images, 
the effect of bad pixels should be negligible. 

• We have ignored the effects of galaxy clus- 
tering, since its amplitude is very weak at these 
magnitudes flNeuschaefer fc Windhorst 1995 , Roche 
et Id. 1993|) . 

• In reality, different fields have different values 
for the seeing, galaxy and star number counts, 
etc., than those used in the simulated data. How- 
ever, the Olw field that we simulated is fairly typ- 
ical in these respects, so generalizing the results 
of the simulated Olw data to the entire survey 
should be valid. 

• We have not taken into account that one of 
our fields, the Ole field, has only one offset im- 



those of the other fields. This will likely raise 
the contamination fraction and possibly reduce 
the completeness slightly, but since the area af- 
fected is <10% of the total the effect should be 
negligible. 

• One possibly significant effect not included is 
the variable PSF over many of our images, due 
mostly to a softening of the focus at increasing 
distances from the center of the field. This will in 
general have the effect of increasing the contam- 
ination, since the classification of faint objects is 
more difficult when the PSF is broader. How- 
ever, only a small fraction of the total area of 
the survey has an extremely broad PSF, so its 
effects should not greatly alter the conclusions of 
our simulations. 

The most important potential problem with 
our simulations is the scale radii used for the 
simulated galaxies in them. Recent results from 
the HST Medium Deep Survey (MDS; [Griffiths 
et al. 1994, Im et al. 1995| ) suggest that many 



agtj in each filter, so thai ils completeness and 
contamination properties will be different than 



faint galaxies have very small half-light radii (me- 
dian 0?4 at V=21-24). Naturally, galaxies with 
small scale radii are more easily misclassified as 
stars than galaxies with large scale radii under 
any seeing conditions. Our simulated galaxies 
have half-light radii from , .'5-l'.'5 at faint mag- 
nitudes; thus, we may have underestimated the 
contamination fraction in the stellar object cat- 
alog. Based on the fraction of our simulated 
galaxies misclassified as stars function of 

their scale radii, we estimate that using the MDS 
scale radii distribution would lead to a doubling 
of the contamination fraction in the stellar ob- 
ject catalog. The contamination fraction is also 
dependent, albeit to a lesser degree, upon the 
relative numbers of galaxies with de Vaucouleurs 
and exponential profiles, and the inclination an- 
gle distribution of the latter population. How- 
ever, these are second-order effects, and our cho- 
sen values for the parameters involved should be 
valid for all but the faintest magnitudes (cf. |im| 



et al. 1995) 



The MDS results are preliminary, being based 
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on a small number of galaxies; also, due to the 
small area studied so far, there is very little 
data on the distribution of galaxy scale radii at 
brighter magnitudes. It is thus perhaps not sur- 
prising that results from the Canada- France Red- 
shift Survey show that only 1% of all galaxies at 
17.5<Iab<22.5 (~19.5<B<24.5) are as compact 



as stars on images with Of! 9 seeing (Crampton 
et ah 1995 ). While our seeing is rarely that 



good, this result indicates that the distribution of 
galaxy scale radii used in the simulations may not 
be as unrealistic as the MDS results suggest, par- 
ticularly at brighter magnitudes. In the next few 
years more extensive data from the MDS should 
yield very accurate distributions of galaxy scale 
radii, morphological types, and axis ratios over a 
large magnitude range. This will enable more ac- 
curate simulations to be done to help understand 
the results of ground-based imaging. Based on 
the preliminary MDS results, we feel that an up- 
wards adjustment of 30% in the contamination 
fraction of the stellar object catalog is reason- 
able. (Recall that this adjustment results in a 
excellent fit to the observed excess of faint stel- 
lar objects above an extrapolation of the counts 
at brighter magnitudes.) 

Thus, in summary, we adopt as the total con- 
tamination fraction of our catalog the prediction 
of our simulated data, scaled to the expected 
number density of faint galaxies and the size 
of our survey, then adjusted upwards by 30%. 
Since the catalog contains 21375 objects, its over- 
all contamination fraction is 40% (8535 objects), 
with roughly half of that occurring fainter than 
the 5a limit. We parametrize the integral con- 
tamination in each filter (i.e. the total number 
of misclassified galaxies in the catalog brighter 
than a given magnitude in that filter) as a cu- 
bic function rising from zero at the threshold 
magnitude to its maximum value at the 3a limit- 
ing magnitude. These magnitudes are tabulated 
for each filter in Table 3, along with the inte- 
gral contamination at the 5a limiting magnitude. 
The completeness of the catalog in each filter 



is 100% down to the threshold magnitude and 
then declines linearly to 90% at the 5a limiting 
magnitude and 50% at the 3a limiting magni- 
tude. The contamination (~25%) and complete- 
ness (90%) of the stellar object catalog at the 5a 
limits are well characterized and allow for reliable 
corrections to be made when computing number 
counts, etc. 

4.4. Photometry and Aperture Correc- 
tions 

Aperture photometry of all stellar objects was 
performed with the IRAF phot task, using a 
centroiding algorithm to determine the position 
of the aperture center. The instrumental magni- 
tudes were calculated using a 3-pixel (1'.'6) radius 
circular aperture and a sky annulus 2'.'7 wide 8" 
away. This was found to be the best compro- 
mise between maximizing the enclosed signal and 
reducing the scatter between different measure- 
ments of the same object due to focus and seeing 
variations between different exposures. Aperture 
photometry was performed at the object's esti- 
mated coordinates in all filters, even if the object 
had not been originally detected in the particu- 
lar filter, in order to obtain color limits on such 
objects. In such cases, occasional centering er- 
rors would occur and need to be corrected - 
the centering routine would latch on to a nearby 
object or noise spike and center the aperture in 
that position, rather than the object's expected 
position. 

Aperture corrections were used to correct the 
magnitude measured in a 1'.'6 aperture to the 5'/8 
aperture used for the standard stars. The correc- 
tions were interactively calculated for the 9 sub- 
sections on each image separately, using up to 20 
of the brightest stars in each subsection. The fi- 
nal variations of the aperture corrections across 
the chip were found to be acceptably small and 
to yield good independent magnitude measure- 
ments from the two offset frames. 
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4.5. Final Magnitude Determination 

Before the determination of final magnitudes, 
the instrumental magnitudes were converted to a 
standard system, as detailed in §5. Once the cal- 
ibrated magnitudes were available for each offset 
image, the final magnitude for each object was 
determined as follows: 

• use the average of the two measurements, un- 
less one measurement lies on a bad pixel, is sat- 
urated, or is a limit 

• if the object lies on a bad pixel or is saturated 
in both images, use the average, but flag it as 
possibly being in error 

• if measurement A is a limit, and B is a detec- 
tion below frame A's limiting magnitude, take 
the detected value from B as the magnitude 

• if measurement A is a limit, and B is a detec- 
tion above frame A's limiting magnitude, average 
the two and set the error negative (this last step 
is done to help mitigate the effects of spurious 
detections caused by overlooked bad pixels) 

• if the final error from the previous steps is 
greater than 0.333, i.e. the object is fainter than 
our 3<7 detection limit, replace the magnitude of 
the object with the 3a magnitude limit of the 
combined frames, and the error with 0.333 

The average 3a magnitude limits in our survey 
are given for each filter in Table 3. Technically, 
3a detections only require a <0.333, regardless 
of the object's magnitude, but it is useful to know 
the average magnitudes to which this limit cor- 
responds. Also, upper magnitude limits for each 
filter, established by examining the magnitudes 
of saturated objects, are given in Table 3. Only 
a handful of saturated objects were found fainter 
than those limits. 

5. Photometric Calibration 

Standard stars in M67, M92, and NGC7790 
were observed on all nights of both observing 
runs. Standard magnitudes and color indices 
were taken variously from Johnson and Sandage 



(1955), Eggen and Sandage (1964), Sandage (1966), 
Schild (1983), Christian et al. (1985), and Ode- 
wahn et al. (1992). These data, supplemented 
with unpublished photometry from A. Porter and 
L. Davis, are available through the IRAF phot- 
CAL on-line library. 

Standard star images were reduced in the nor- 
mal manner for CCD images, using the same 
domeflats and skyflats used for the field images. 
Instrumental magnitudes were found for all stan- 
dard stars using 5'.'8 radius circular apertures, 
with all other parameters, including the sky an- 
nulus size, being the same as used for the field 
images. 

The photometric solution for the UBVR filters 
is independent of that for our two I filters, and 
is discussed first. In the former case we trans- 
form our observed magnitudes to the standard 
UBVR system; in the latter we define a new sys- 
tem for the I filters. 

5.1. UBVR Calibration 

The IRAF photcal photometric calibration 
package was used to interactively solve the fol- 
lowing system of equations involving the instru- 
mental magnitudes (lower case), the standard 
calibrated magnitudes and colors (upper case), 
and the airmass for the observations in each fil- 
ter (X u , Xf,, X„, X r ). The system of equations 
was solved simultaneously for the zeropoints (uo, 
bo, vo, ro), the color coefficients (ui, bi, vi, ri), 
and the extinction coefficients (u2, b2, v 2 , r2) for 
each night separately. 

u = u +u 1 x(U-B)+u 2 xX u +V+{B-V)+{U-B) 

(1) 

b = b + b 1 x(B-V) + b 2 xX b + V + (B-V) (2) 
v = v + vi x (V-R) + v 2 x X v + V (3) 
r = r + r 1 x(V-R) + r 2 xX r + V-(V-R) (4) 

Not all nights yielded well-defined photomet- 
ric solutions independently, since typically only 
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three observations of standard star clusters were 
made each night. Since most parameters in the 
solutions had the same values within the errors 
from night to night, and indeed between the 
two runs, we decided to adopt the same color 
terms and extinction coefficients for each observ- 
ing run and let only the zeropoint in each filter 
vary. Thus we solved the equations again, this 
time only for the zeropoint, keeping the other 
variables fixed at the average values determined 
from the individual nights' data. The resulting 
final UBVR photometric calibration parameters 
for both runs are listed in Table 4. 

As shown in Table 4, the formal errors associ- 
ated with the photometric parameters £1X6 8l few 
percent at most. However, differences from the 
base zeropoints were found to exist for most CCD 
images of the survey fields in both runs, as dis- 
cussed in §5.3. 

5.2. 175 and 186 Calibration 

Our magnitude system in the 175 and 186 fil- 
ters is defined by observations of the spectropho- 
tometry standard DA white dwarf Wolf 1346 
made on 1990 August 24. The definition and 
calibration of our system to the AB system are 
described in detail in this section. 

Instrumental magnitudes for Wolf 1346 were 
determined using a (8'/5) radius circular aper- 
ture. AB magnitudes for the star were taken 
from the IRAF onedstds on-line library, which 
in the case of Wolf 1346 consists of unpublished 
data plus data from Oke (1974) for wavelengths 
redward of 8400 A. Despite its greater noise, we 
chose to use the Oke data rather than the data 
from Massey et al. (1988) because of its better 
coverage of the 7000-9500 A wavelength region. 

In order to find the effective wavelengths of 
the filters for this star, the AB magnitudes were 
converted to fluxes and were convolved with the 
filter + CCD + telescope system response curves. 
The effective wavelengths for Wolf 1346 in the 
175 and 186 filters are 7430 A and 8580 land 



the standard magnitudes of Wolf 1346 are thus 
defined to be its interpolated AB magnitudes at 
these wavelengths: 175=12.02 and 186=12.22. 

The extinction coefficients at the wavelengths 
of these filters was calculated in two ways: first 
by using the standard KPNO extinction curve 
from the IRAF onedstds library to find the 'ef- 
fective extinction coefficient' and second by com- 
paring measurements of the instrumental magni- 
tudes of stars in NGC 7790 made at different air- 
masses during the same night of the 1990 August 
run. The extinction coefficients for 175 found by 
the standard and empirical methods were 0.063 
and 0.078±0.014 respectively; and for 186 they 
were 0.039 and 0.039±0.008. We chose to adopt 
the empirical values, but the difference between 
using the different values would be only a few 
percent at most, since 94% of our observations 
were made at airmasses less than 1.4 and 100% 
at airmasses less than 1.65. 

The final step in determining the zeropoints is 
simply one of plugging in to the equation 

fninstr = fnAB + zeropoint + X * airmass (5) 

which defines the zeropoints of the 175 and 186 
filters for the 1990 August run. These zeropoints 
were found to be -22.08 and -21.53, as given in 
Table 4. Note that no color term was adopted 
for our two I filters. 

In order to obtain the 175 and 186 filter zero- 
points for the 1991 April run, the standard stars 
in the standard fields observed during both runs 
(M92 and NGC7790) were defined as secondary 
standards in the 175 and 186 filters as follows. 
Using the data from 1990 August, equation (5) 
was solved for uiab for each standard star in 
each cluster, holding fixed the zeropoints and the 
extinction coefficients at the values given above. 
The values of thab derived from these 1990 Au- 
gust observations were adopted as the magni- 
tudes for these standard stars in our two I fil- 
ters. The zeropoints for the 1991 April runs were 
then found by solving the standard photometric 
transformation equations in the same manner as 
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equations (1-4) for the UBVR niters, with the 
exception that no color terms were used. Once 
again, the extinction coefficients were held fixed, 
because we assumed that the extinction did not 
change markedly between the two runs, which 
was certainly the case with the UBVR data. 

i75 = 775 + i75 + *75 2 x X i75 (6) 
i86 = 786 + i86 + i86 2 x X m (7) 

These zeropoints were found to be -22.03 and 
-21.48, as given in Table 4. These values are only 
5% different from those for 1990 August. Again, 
as for the UBVR filters, some differences from the 
base zeropoints in were found to exist for most 
CCD images of the survey fields in both the 1990 
August and 1991 April runs, as listed in Table 2. 

5.3. Further Empirical Calibration 

The photometric calibration derived from the 
standard stars was checked on the survey fields 
by calibrating the instrumental magnitudes us- 
ing only the zeropoint and airmass coefficient at 
first, i.e. without including the color term. This 
allowed us to check for any zeropoint shifts be- 
tween CCD images of the same field and make 
adjustments to the zeropoints as follows. We 
plotted the difference between the magnitudes 
measured for all 'good' objects (ones not af- 
fected by bad pixels) in the two independent 
offset frames of each field versus the average of 
the two magnitudes (hereafter we refer to these 
as the 'magnitude difference' graphs). A typi- 
cal magnitude difference graph is shown in Fig- 
ure 3. Figure 3 shows that at bright magnitudes 
the two measurements agree to within the small 
(<2%) uncertainty due to limitations of the flat- 
fielding, etc. It also shows that our bad pixel 
removal is quite good, with only 2 or 3 objects 
(out of ~1000) showing differences between the 
two magnitude measurements large enough to be 
undisputably due to errors. If the average mag- 
nitude difference in a plot like Figure 3 is not 
zero, one of the frames has been calibrated with 
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Fig. 3. — A typical magnitude difference graph, 
in this case from the V images in the Olw field. 
The x-axis is the average V magnitude measured 
in the two images for each object, and the y-axis 
is the difference between the magnitudes mea- 
sured in the two frames for all 'good' objects 
(ones not affected by bad pixels or focus prob- 
lems) . 

an erroneous zeropoint. The magnitude differ- 
ence graphs are only useful in calculating zero- 
point shifts when the same set of objects has 
been measured; so, in order to measure zeropoint 
shifts between frames in different fields, we plot- 
ted and compared various color-color diagrams 
as detailed below. 

It became apparent from these graphs that 
empirical corrections to the photometric solu- 
tions in most of our fields were necessary in order 
to bring all fields onto a common system, similar 
to the practice in astrophotography of calibrating 
each photographic plate independently. There 
are many possible explanations for zeropoint dif- 
ferences between fields - the possible presence 
of light cirrus as noted in our observing logs, 
as well as night-to-night variations in the atmo- 
spheric transparency and extinction, being the 
most prominent - but for our purposes the exact 
causes are not relevant. This is because one of 
the main reasons we calibrate the photometry is 
to bring all the fields to a common system where 
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the locus of stellar colors shows as little scatter as 
possible, which helps us better discriminate can- 
didate quasars from stars. It is not as efficient 
to search for outliers in each field independently 
because of the relatively small number of objects 
that define the stellar locus in each field. That 
is, a slightly unusual but normal star (e.g. O, B, 
and A stars at high galactic latitude) will stand 
out more in the color-color diagrams of each in- 
dividual field than in the combined diagram of 
all fields. Strictly speaking, our empirical correc- 
tions should be made to the extinction terms in 
the photometric solution, since they are the most 
likely to vary. However, it was easier to derive a 
correction to the zeropoint of each field rather 
than to the extinction term, so this was done in- 
stead. Since the shift needed to align the stellar 
locus of a given field is fixed, the term used to 
make that adjustment matters little in practical 
terms. 

The 17n field was adopted as our best field, 
relative to which all others were adjusted. It 
showed an average magnitude difference between 
offset frames of zero in all but the U and B filters, 
and in those filters it was obvious by inspection 
and by consulting our observing logs which two 
frames were affected by light cirrus. Thus we 
believe that the 17n magnitude calibration is ex- 
act to within a few percent. Note, however, that 
since we have tied our magnitude calibration to 
this field, any systematic error in the calibration 
assumed for it will affect all our other fields as 
well. Since the 17n field was observed in the 
1991 April run, all other 1991 April fields had 
their zeropoints adjusted relative to the 17n field 
by direct comparison of several color-color dia- 
grams. For both sets of offset frames in each 
field, color-color diagrams were plotted on an 
identical scale using the brightest one hundred 
unsaturated non-bad-pixel objects. Each field's 
diagrams were compared with those for the 17n 
field and were shifted by eye in both colors until 
the best alignment was obtained. The empiri- 
cally required shift in a given color is simply a 



combination of the empirical zeropoint shifts re- 
quired for the two filters comprising the color. 
By combining the shifts required for five (inde- 
pendent) colors with the shifts required between 
frames of the same filter in each field, we are able 
to determine each field's zeropoint shifts relative 
to the 17n field to within an additive constant 
which is the same for all filters. That is, we can 
determine the relative zeropoint shifts between 
all filters, but there is still an overall additive 
constant for the zeropoint in all filters that we 
cannot determine in the same manner. This ad- 
ditive constant was determined by requiring that 
all empirical R-band zeropoint shifts be as small 
as possible relative to the 17n field. This is prob- 
ably a better estimate than requiring the average 
zeropoint shift in all filters to be minimized, be- 
cause the U, B, and V filters will be affected by 
varying atmospheric extinction more than the R 
band. For our 1990 August run, we used the fact 
that the 17s field, which was observed at that 
time, overlaps with the 17n field. Approximately 
twenty of the brightest (R <20) objects in the 
overlap region were used to find the exact shifts 
between the zeropoints for the two fields by re- 
quiring that the average difference between these 
two independent measures of the magnitudes of 
these objects be equal to zero. The zeropoints of 
17s were adjusted as needed to achieve this, and 
then the other 1990 August fields were adjusted 
relative to the 17s field in the same manner as the 
1991 April fields were adjusted relative to the 17n 
field. The average shifts in each filter were less 
than 5% half the time, and less than 10% in all 
cases except the R filter in the 1990 August run, 
where the empirically determined shift was 0.13 
magnitudes, or 13%, different from its value cal- 
culated by the photometric comparison. While 
this is uncomfortably high, comparison of color- 
color plots from different fields shows clearly that 
our empirical zeropoint shifts achieve the desired 
result of tightening the stellar locus. Since the 
empirical shifts in the zeropoints are typically be- 
tween 5 and 10%, we estimate that the system- 
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atic uncertainties in our magnitude calibration 
are at this level. 

As a final check, we used bright (R< 20) ob- 
jects, typically two dozen of them, in the overlap 
regions between the 21-hour and 22-hour fields 
to confirm our determinations of the zeropoint 
shifts for the different filters and fields. The mag- 
nitudes agreed to within 2a in all cases. The fi- 
nal zeropoints for each CCD field in both runs 
are given in Table 2. Note that for all filters in 
the Ole field and the B filter in the lOe field, only 
1 image was obtained. Thus no magnitude dif- 
ference plot could be constructed, and calibra- 
tion of the zero points was performed only by 
color-color diagram comparison. This should not 
introduce any significant additional uncertainty 
into the calibration of these fields. 

5.4. Photometric Error Analysis And Dis- 
cussion 

We also used the magnitude difference graphs 
to study the accuracy of the photometric errors 
determined by the IRAF photometry routines, to 
ensure that there were no unexpected additional 
sources of error that had not been taken into ac- 
count. 

The magnitude difference is d= mi-m2, where 
mi±(Ji and m2±o"2 are values of the magnitude 
of the same object measured from independent 
CCD frames. Thus the distribution of magnitude 
differences at any particular magnitude should be 



given by a gaussian: 

f(d(m avg )) = (l/fj d v / 27r)e 



-[d(m)] 2 /2^ 



(8) 



where the magnitude m a „ 9 =(mi+ni2)/2 and the 
sigma of this distribution is (T^(m avg ) = a\ + o\. 
It is reasonable in most cases to assume o\ = oi , 
i.e. that the photometric errors in each CCD 
frame of a particular filter and field are the same 
function of magnitude. Thus if we measure 
0d, we can find a from the equation a = a d j 

We obtained two different estimates of <Jd from 
our magnitude difference plots. First, we took 



the square of the magnitude difference for every 
object not affected by a bad pixel, sorted by mag- 
nitude, binned them in groups of 50, and calcu- 
lated the average value of the square of the mag- 
nitude difference for each bin, Second, we took 
the absolute value of the magnitude difference 
for each object, sorted and binned them as be- 
fore, and calculated the average absolute value in 
each bin. This gives us a measure of the mean 
deviation a d = J \d\f(d)=0.8a d flYoung 1962]) . 



In Figure 4 we plot the photometric errors 
from the IRAF phot routine for all objects in 
both CCD frames in a representative field and 
filter and the sigma estimated from the magni- 
tude difference plots in both manners described 
above. The match is clearly quite good. At the 
brightest magnitudes we have replaced the pho- 
tometric errors determined by IRAF with 0.02 if 
they were originally smaller than this (to be ex- 
act, this error should be added in quadrature to 
all the objects' errors). 

In summary, the external errors in our pho- 
tometry, as estimated from the two independent 
measurements of the magnitudes of each of our 
objects, have been shown to be consistent with 
the internal errors computed according to pho- 
ton statistics, except for a ~2% additional uncer- 
tainty independent of magnitude. This indepen- 
dent check proves that our flatfielding process, 
aperture correction procedures, and photometry 
methods are all quite reliable, having inherent 
limitations of only the aforementioned ~2%. As 
for systematic errors, our stellar locus matches 
values for stellar colors from the literature to 
about 5%. An exact match with expected col- 
ors is unimportant in terms of outlier selection, 
but necessary for accurate comparison with other 
studies and for other uses of the dataset. 

6. Astrometry 

The final step in the preparation of the cata- 
log was the conversion of CCD (x,y) coordinates 
(centroids) to right ascension and declination. 
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We required accuracy to <0 / /5 because the op- 
tical fibers used by the HYDRA instrument for 
our initial spectroscopy are only 2'.'0 in diameter. 
The astrometry was performed by A. Porter in 
the manner briefly outlined here. 

The program FINDER ( Anderson 199p was 
run on each survey field to search for HST Guide 
Star Catalog (GSC; |Taff et al. 1990|) stars within 
the field. There were sometimes less than a dozen 
or so GSC stars in a field. In these cases we 
used POSS plates to construct an astrometric 
grid based on nearby AGK stars. The positions 
of all the GSC stars in the field, plus a dozen or 
two anonymous field stars, were then measured 
on this grid. This served check on the GSC 
astrometry (which was found to be accurate to 
<0f!b in our fields, adequate for our purposes), 
and yielded enough stars (between 14-26 for each 
CCD field) to construct a good astrometric solu- 
tion. The CCD (x,y) centroids for all these as- 
trometric calibration stars were measured on the 
first R frames of each field using the IRAF task 
imexamine. The astrometric solutions for CCD 
(x,y) to (RA,dec) conversion in each field were 
calculated using the program ASTRO flAnder 



sor 1992) and included terms up through a r 
barrel distortion term. The internal residuals of 
the solutions ranged from (/.'18 to (//46, and no 
significant trends or distortions in the solutions 
were present. The average positional agreement 
between stars in the overlap regions of adjacent 
CCD fields was usually better than (/.'l and al- 
ways better than Of! A. The solutions were then 
used to find the RA and Dec of all objects in the 
catalog, using the CCD (x,y) positions of the ob- 
jects as measured from the first R frame of each 
field. 



errors are consistent with expectations from pho- 
ton statistics. The production of the catalog is 
the main result of this paper. From this catalog 
we would like to isolate as many as possible of the 
~193 quasars with B<22.6 expected in our sur- 
vey area (Majewski et al. 1993), with a minimum 
of contamination from stars and other objects. 
As there are 4,205 objects with B<22.6, quasars 
are outnumbered by a factor of 20:1, and careful 
candidate selection is obviously necessary. Paper 
II flHall et al. 1996; ) discusses our candidate selec- 
tion methods and presents initial spectroscopic 
results and conclusions. 
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ticularly the telescope operators, for their assis- 
tance in observing, and last but definitely not 
least the KPNO TAC for their allocation of time 
for this project. PBH acknowledges support from 
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7. Summary 

The final photometrically calibrated catalog 
of stellar objects in our Deep Multicolor Survey 
contains 21,375 objects from 0.83D of sky. All 
objects have positions accurate to ~0 / /5 in each 
coordinate and magnitudes in six bands whose 
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Table 1 

Deep Multicolor Survey Fields 



Field 


RA (1950.0) 


Dec (1950.0) 


1(H) 


b(II) 


Area a 


E(B-V) b 


Ole 


01:00:16 


-00:57:52 


129.3069 


-63.4266 


292.4 


0.0192 


Olw 


00:59:21 


-00:57:55 


128.7997 


-63.4488 


224.9 


0.0273 


lOe 


10:34:19 


-00:44:40 


248.4348 


47.1235 


286.0 


0.0302 


14n 


13:58:40 


-00:37:06 


336.7107 


57.3613 


288.5 


0.0283 


14s 


13:58:40 


-00:54:24 


336.4456 


57.1110 


272.8 


0.0263 


17n 


17:14:43 


+50:16:47 


76.9826 


35.5101 


266.0 


0.0082 


17s 


17:14:45 


+50:00:01 


76.6398 


35.4963 


285.9 


0.0093 


21e 


21:39:57 


-04:00:06 


51.8119 


-39.2835 


247.0 


0.0243 


21w 


21:39:01 


-04:00:01 


51.6505 


-39.0849 


290.0 


0.0203 


22c 


22:47:04 


-02:09:11 


68.6033 


-51.4310 


286.9 


0.0542 


22w 


22:46:07 


-02:09:14 


68.3400 


-51.2595 


249.7 


0.0492 



a Units of square arcminutes. 

b Calculated with a program supplied by D. Burstein, using data given in Burstein & Heiles 1978 and Burstein 
& Heiles 1982. 



Table 4 

Photometric Calibration Parameters 





Color 


Extinction 


1991 Base 


1991 Average 


1990 Base 


1990 Average 


Filter 


Term 


Coefficient 


Zeropoint 


Zeropoint 


Zeropoint 


Zeropoint 


U 


-0.043+0.011 


0.530+0.023 


-20.64 


-20.69 


-20.78 


-20.71 


B 


-0.063+0.016 


0.275+0.015 


-21.71 


-21.75 


-21.80 


-21.74 


V 


-0.092+0.018 


0.113+0.007 


-22.45 


-22.50 


-22.51 


-22.46 


R 


-0.098+0.020 


0.087+0.017 


-22.06 


-22.06 


-22.16 


-22.03 


175 




0.078+0.014 


-22.03 


-22.05 


-22.08 


-21.99 


186 




0.039+0.008 


-21.48 


-21.48 


-21.53 


-21.45 
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Table 2 
Observations 



Field Filter 


Frame 


Exposure a 


UT Date 


UT Time 


Airmass 


Seemg( ) 


Zeropoint 


nip tt 

u±e \j 




1 


on /8 /94 
yu /of 


1 1 -1 % 
ii.it> 


1 91 Q 


1 RO 


90 RQ 
-zu.uy 


B 


1 


1080 


90/8/93 
yu I o / Zj'j 


11:00 


1.194 


1.21 


-21.66 


v 


\ 


^40 


qn /ft/9^ 
yu / of £Q 


1 0-4^ 


1 1 77 

1.11 1 


1 10 

1.1U 


99 ^7 


R 


\ 


^40 


CfU j Of ZrO 


11:22 


1.216 


1.22 


-22.02 


1 i o 


1 


^40 


on /S./9*} 
yu / of ZrO 


1 1 

1 1 .00 


1 9^0 


1.28 


91 Q8 
ii. yo 




i 

1. 


i 080 

1UOU 


on lR /OA 


1 1 -41 
ii.ii 


1 9fi8 


9 0^ 


91 48 




\ 


i ^o 

louu 


on /ft/99 


1U. 1J 


1 1 70 

1.1IU 


1 1 


90 R<\ 


TT 




i ^o 


qn/8/99 
yu /Of 


10:41 


1 1 78 
1 . 1 1 


1.24 


-20.82 


B 


x 


1080 


Q0/8/9S 
yu /of Zj'j 


10:09 


1.170 


1.53 


-21.77 


B 




i oso 

1UOU 


qn/8/91 
yu / o / £±. 


10:42 


1 1 7^ 
1 . 1 1 


1 

i .uy 


99 7Q 


v 


x 


540 


QD/8/9S 
yu o 


10:31 


1.172 


1.26 


-22.46 


v 


o 
z . 


^40 

OIU 


qn /o /oi 
yu/ o/ Zr± 


1 0-9Q 


1 1 70 
J — i / u 


1 1 4 

1.11 


99 47 


R 


\ 


540 


QO/8/99 
yu o 


9:48 


1.181 


0.93 


-22.02 


1 1 


2. 


^40 

UlU 


qn/s /91 
yu j of & -L 


11:04 


1 1 87 
1.101 


1 OR 

l.UU 


91 0" 1 ! 


1 / u 


i 

i . 


^40 
oiu 


qn /o /oi 
yu / o/ Zr± 


1 0-09 


1 17Q 
1.1/0 


1 00 

l.UU 


91 Qfi 
-z i . yu 


175 


2. 


540 


QD/8/91 
yu f o f ^ i. 


11:18 


1.201 


0.99 


-21.99 


186 




1080 


90/8/99 
yu f o j 


11:31 


1.228 


1.39 


-21.44 


TSfi 


2, 


1 080 
1UOU 


qn/8/99 

yU / O/ ZiZr 


11:14 


1 907 




u. yu 


-21.41 


lOe U 




1350 


91/4/15 


6:07 


1.338 


1.36 


-20.74 


TT 


2, 


1 

louu 


Q1 /4/1 ^ 

y _1_/ 4:/ LO 


fi-4^ 


1 47 c > 


1 ^7 


90 7^ 


B 


1. 


1080 


91/4/15 


7:10 


1.615 


1.29 


-21.73 


V 


1. 


540 


91/4/14 


4:34 


1.191 


1.41 


-22.49 


V 


2. 


540 


91/4/14 


4:47 


1.195 


1.62 


-22.48 


R 


1. 


675 


91/4/14 


5:02 


1.204 


1.57 


-22.06 


R 


2. 


675 


91/4/14 


5:17 


1.218 


1.66 


-22.06 


175 


1. 


540 


91/4/14 


5:32 


1.237 


1.69 


-22.05 


175 


2. 


540 


91/4/14 


5:45 


1.258 


1.72 


-22.05 


186 


1. 


1080 


91/4/14 


6:00 


1.302 


1.65 


-21.48 


186 


2. 


1080 


91/4/15 


5:43 


1.244 


1.02 


-21.50 


14n U 


1. 


1350 


91/4/14 


9:08 


1.269 


2.35 


-20.71 


U 


2. 


1350 


91/4/14 


9:35 


1.336 


2.66 


-20.77 


B 


1. 


1080 


91/4/15 


7:45 


1.190 


1.10 


-21.77 


B 


2. 


1080 


91/4/15 


8:07 


1.195 


1.16 


-21.76 


V 


1. 


540 


91/4/14 


7:12 


1.210 


1.54 


-22.50 


V 


2. 


540 


91/4/14 


6:59 


1.225 


1.48 


-22.50 


R 


1. 


675 


91/4/14 


7:26 


1.198 


1.74 


-22.06 


R 


2. 


675 


91/4/14 


7:41 


1.191 


1.70 


-22.06 


175 


1. 


540 


91/4/14 


8:09 


1.192 


1.95 


-22.06 


175 


2. 


540 


91/4/14 


7:56 


1.189 


2.03 


-21.06 
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Table 2 — Continued 



Field Filter 


Frame 


Exposure a 


UT Date 


UT Time 


Airmass 


Seemg( ) 


Zeropoint 


TSfi 
±ou 


1 


1 080 

1UOU 


Q1 /4/14 


8-9^ 
o.zo 


1 909 
i .zuz 


1 87 


91 48 


TSfi 


2, 


1 080 

1UOU 


Q1 /4/14 
y i / 1 / n 


8-4^ 


1.225 


1 Q1 

i .y i 


-21.48 


1 4d TT 


1 


1 Vi,C\ 


Q1 /4/1fi 


y.jy 


1 4 c iQ 
i .uy 


1 fiQ 


90 fil 

ZU.U1 


TT 

U 


o 
z . 


1 9^0 


Q1 /4 /I fi 


8-4^ 


1 944 
i .zii 


1 4** 

1 .10 


90 ^7 


T3 
JZ> 


1 

i . 


1 080 

1UOU 


Q1 /4 /I 


7-10 


1 908 
1 .zuo 


1 07 


91 77 

-z 1 . t 1 


[ J 


o 
z . 


1 080 

±UOU 


Q1 /4 /I 7 


Q-48 


1 490 


1 88 

L.OO 


91 77 

-Z 1 . ( J 


v 

V 




^40 


Q1 /4/1fi 
y 1/ *±/ iu 


8-01 
o.ui 


1 1 Qfi 
i . lyu 


1 ^1 


-22.30 


V 

V 




^40 


Q1 /4/1 'i 

y 1/ 4:/ 1U 


O.UI 


1 9^8 

1 .ZOO 


1.02 


-22.50 


R 
XX 


i 

_L . 


fi7<1 


Q1 /4/1 fi 
y i/ 1/ iu 


O.Ol 


1 .zuo 


1 40 

1 .IU 


91 33 

-Zi L.OO 


R 


9 
z. 


fi7<i 


Q1 /4 /I ^ 


Q-O^ 
y.uo 


1 9fi ,: t 
1 .zuo 


1 08 


99 Ofi 
-zz. uu 


1 ( O 




^40 


01 /4/1 7 


O.O-L 


1.221 


2.12 


-22.03 


L 1 o 


2, 


^40 


Q1 /4/1fi 

y 1/ 1/ IU 


S-^l 


1.216 


1.44 


-21.71 


TSfi 




1 080 

1UOU 


Q1 /4/14 


fi-00 
u.uu 


1 9Q7 
i .zy i 


1 1 7 


-21.46 


TSfi 


2, 


1 080 

1UOU 


Q1 /4/1fi 
y 1/ 1/ iu 


y.OU 


1 9 c iQ 


1 97 

1 .Z ( 


-21.44 


17n U 


x 


1350 


91/4/16 


10:47 


1.053 


1.57 


-20.38 


TT 


2, 


1 "^0 
louu 


Q1 /4/1 ^ 
y i / 1/ 10 


11:01 


1 O 1 ^ 


1 Ofi 

l.UU 


90 P,<\ 
zu.uu 


J3 




1 080 

1UOU 


Q1 /4/1 'i 
y 1/ 1/ ±o 


1 0-H7 

1U.O 1 


1 OW 


Q7 


91 73 

Z 1 . 1 o 


B 


2. 


1080 


91/4/16 


11:14 


1.055 


1.40 


-21.52 


V 

V 


\ 


^40 


Q1 /4/14 
y i / 1 / n 


10:21 


1 Ofifi 
J- .uuu 


1 47 


-22.50 


V 

V 


2, 


^40 


Q1 /4/14 

y 1 / 1 / 11 


1 0-^fi 

1U.OU 


1 O^Q 


1 47 


-22.50 


R 


1. 


675 


91/4/14 


10:50 


1.055 


1.34 


-22.06 


R 


2. 


675 


91/4/14 


11:06 


1.053 


1.25 


-22.06 


175 


1. 


540 


91/4/14 


11:21 


1.054 


1.20 


-22.06 


175 


2. 


540 


91/4/14 


11:33 


1.056 


1.16 


-22.06 


186 


1. 


1080 


91/4/15 


11:29 


1.058 


0.90 


-21.48 


186 


2. 


1080 


91/4/15 


10:14 


1.065 


0.90 


-21.48 


17s U 


1. 


1350 


90/8/21 


5:42 


1.289 


1.27 


-20.76 


U 


2. 


1350 


90/8/24 


4:20 


1.137 


1.26 


-20.75 


B 


1. 


1080 


90/8/21 


4:53 


1.169 


1.39 


-21.79 


B 


2. 


1080 


90/8/23 


4:04 


1.107 


1.44 


-21.77 


V 


1. 


540 


90/8/21 


4:40 


1.139 


1.11 


-22.49 


V 


2. 


540 


90/8/23 


4:26 


1.129 


1.15 


-22.49 


R 


1. 


540 


90/8/23 


4:45 


1.159 


1.21 


-22.02 


R 


2. 


540 


90/8/21 


5:16 


1.205 


1.23 


-22.05 


175 


1. 


540 


90/8/21 


5:29 


1.234 


1.06 


-21.99 


175 


2. 


540 


90/8/23 


4:57 


1.183 


1.35 


-21.95 


186 


1. 


1080 


90/8/23 


3:41 


1.083 


1.25 


-21.40 


186 


2. 


1080 


90/8/23 


6:10 


1.371 


1.30 


-21.41 
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Table 2 — Continued 



Field Filter Frame Exposure a UT Date UT Time Airmass Seeing(") Zeropoint b 



u 


1. 


1350 


90/8/24 


u 


2. 


1350 


90/8/24 


B 


1. 


1080 


90/8/23 


B 


2. 


1080 


90/8/24 


V 


1. 


540 


90/8/24 


V 


2. 


540 


90/8/23 


R 


1. 


540 


90/8/24 


R 


2. 


540 


90/8/23 


175 


1. 


540 


90/8/24 


175 


2. 


540 


90/8/23 


186 


1. 


1080 


90/8/23 


186 


2. 


1080 


90/8/24 


U 


1. 


1350 


90/8/21 


U 


2. 


1350 


90/8/22 


B 


1. 


1080 


90/8/21 


B 


2. 


1080 


90/8/22 


V 


1. 


540 


90/8/21 


V 


2. 


540 


90/8/22 


R 


1. 


540 


90/8/22 


R 


2. 


540 


90/8/21 


175 


1. 


540 


90/8/21 


175 


2. 


540 


90/8/22 


186 


1. 


1080 


90/8/21 


186 


2. 


1080 


90/8/22 


U 


1. 


1350 


90/8/24 


u 


2. 


1350 


90/8/24 


B 


1. 


1080 


90/8/24 


B 


2. 


1080 


90/8/23 


V 


1. 


540 


90/8/23 


V 


2. 


540 


90/8/24 


R 


1. 


540 


90/8/23 


R 


2. 


540 


90/8/24 


175 


1. 


540 


90/8/23 


175 


2. 


540 


90/8/24 


186 


1. 


1080 


90/8/24 


186 


2. 


1080 


90/8/23 



5:27 


1.318 


1.71 


-20.83 


5:53 


1.271 


1.48 


-20.83 


5:57 


1.274 


1.25 


-21.81 


6:44 


1.233 


1.78 


-21.81 


7:06 


1.234 


1.54 


-22.50 


5:42 


1.308 


1.28 


-22.49 


7:19 


1.240 


1.57 


-22.03 


6:20 


1.251 


1.21 


-22.02 


7:32 


1.250 


1.38 


-22.00 


6:33 


1.241 


1.34 


-21.99 


7:13 


1.237 


1.11 


-21.44 


6:21 


1.243 


1.28 


-21.44 


7:41 


1.256 


1.86 


-20.66 


6:23 


1.246 


1.44 


-20.72 


6:51 


1.234 


1.86 


-21.77 


5:34 


1.321 


1.19 


-21.77 


6:38 


1.242 


1.55 


-22.49 


5:21 


1.366 


1.17 


-22.48 


5:57 


1.285 


1.02 


-22.02 


7:15 


1.234 


1.56 


-22.06 


7:28 


1.239 


1.52 


-22.01 


6:09 


1.266 


1.00 


-21.99 


6:08 


1.291 


1.64 


-21.48 


6:50 


1.233 


1.54 


-21.46 


8:27 


1.218 


1.80 


-20.62 


8:54 


1.248 


1.38 


-20.66 


9:42 


1.349 


2.32 


-21.68 


8:30 


1.215 


1.06 


-21.69 


8:15 


1.206 


1.09 


-22.44 


10:04 


1.407 


1.86 


-22.45 


8:51 


1.230 


0.98 


-22.07 


10:17 


1.461 


2.05 


-22.00 


9:04 


1.246 


1.11 


-22.01 


10:30 


1.522 


1.83 


-22.03 


9:21 


1.293 


1.63 


-21.48 


9:44 


1.343 


1.09 


-21.46 
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Table 2 — Continued 



Field 


Filter 


Frame 


Exposure a 


UT Date 


UT Time 


Airmass 


Seeing(") 


Zeropoint b 


22w 


U 


1. 


1350 


90/8/24 


7:59 


1.206 


1.57 


-20.52 




u 


2. 


1350 


90/8/22 


9:06 


1.257 


1.23 


-20.46 




B 


1. 


1080 


90/8/22 


8:18 


1.207 


1.24 


-21.62 




B 


2. 


1080 


90/8/21 


9:09 


1.250 


1.59 


-21.67 




V 


1. 


540 


90/8/21 


8:55 


1.227 


1.52 


-22.45 




V 


2. 


540 


90/8/22 


8:05 


1.205 


1.12 


-22.42 




R 


1. 


540 


90/8/22 


8:40 


1.217 


1.01 


-22.02 




R 


2. 


540 


90/8/21 


9:40 


1.302 


1.45 


-22.04 




175 


1. 


540 


90/8/22 


8:53 


1.229 


1.07 


-21.98 




175 


2. 


540 


90/8/21 


9:53 


1.334 


1.54 


-22.00 




186 


1. 


1080 


90/8/23 


7:47 


1.207 


1.11 


-21.47 




186 


2. 


1080 


90/8/22 


7:14 


1.234 


1.57 


-21.47 



a Units of seconds. 
b See §5.3. 



Table 3 
Average Magnitude Limits 





Threshold 


5cr Limiting 


3cr Limiting 


Upper Limit 




Filter 


Magnitude 


Magnitude 8, 


Magnitude 13 


Magnitude 


G(m 5CT ) c 


U 


21.0 


22.8 


23.45 


16.0 


25.2% 


B 


21.5 


23.8 


24.30 


16.5 


22.2% 


V 


21.0 


23.5 


24.00 


16.5 


27.8% 


R 


20.0 


23.0 


23.50 


16.5 


33.0% 


175 


20.0 


22.4 


23.05 


16.0 


19.5% 


186 


19.5 


22.1 


22.70 


16.0 


21.5% 



a Also 90% completeness magnitude. 
b Also 50% completeness magnitude. 

The percentage of misclassificd galaxies in the stellar object catalog at the 5er limiting magnitude. 
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Fig. 4. — Formal photometric errors vs. errors 
estimated from the magnitude difference plots 
like that shown in Figure 3. The points are the 
photometric errors from the IRAF 'phot' routine 
for all objects in both CCD frames in the 21w 
field, 18600 filter. The dashed line is the sigma 
estimated from the average square of the mag- 
nitude difference for every 50 objects, while the 
connected line is the sigma estimated for every 
50 objects using the average absolute value of 
the magnitude difference. The match is clearly 
quite good at all magnitudes. 
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